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Abstract 

A numerical  methodology  for  the  evaluation  of  aero- 
dynamic loads  acting  on  a complex  lifting  config- 
uration is  presented.  The  work  is  limited  to  the 
case  of  attached  high-Reynolds  number  flows . A vis- 
cous/potential interaction  technique  is  utilized  to  take 
into  account  the  effects  of  the  viscosity.  For  the 
potential-flow  analysis,  a boundary  element  formula- 
tion is  used;  for  simplicity,  only  incompressible  flows 
are  examined.  The  theoretical  basis  of  the  present 
methodology  is  briefly  described.  Comparisons  with 
available,  numerical  and  experimental  results  are  in- 
cluded. 

1.  Introduction 

The  present  work  is  an  overview  paper  on  recent  re- 
sults on  viscous-flow  analysis  of  innovative  configu- 
rations (low-induced  drag  lifting  configurations),  ob- 
tained by  the  authors  within  the  objective  of  de- 
veloping an  MDO  methodology.  The  aerodynamic 
methodology  is  based  on  a boundary  element  formu- 
lation for  the  velocity  potential  introduced  by  Morino 
[16];  the  effects  of  viscosity  are  taken  into  account  by 
a classical  viscous/in  viscid  coupling  technique;  a tech- 
nique for  extending  this  approach  beyond  the  limits 
of  the  boundary-layer  approach  is  also  indicated. 

In  order  to  put  this  work  in  the  paper  context,  note 
that,  the  increased  request  for  low-cost  air  trans- 
portation has  generated  considerable  interest  for  non- 
conventional  configurations,  such  as  biplanes  with 
wing  tips  connected  to  each  other  directly  (joined 
wing  configuration)  or  through  a vertical  surface 
(box-wing  configuration).  These  configurations  are 
related  to  the  early  work  of  Prandtl  [25]  on  multi- 


planes, which  shows  that  these  configurations  have 
certain  aerodynamic  advantages  (outlined  in  the  fol- 
lowing sentences)  with  respect  to  isolated  wings. 
First,  the  induced  drag  of  a multiplane  with  elliptic 
distribution  of  circulation  is  lower  than  the  induced 
drag  of  the  equivalent  monoplane  (Le.,  of  the  mono- 
plane having  the  same  span  and  generating  the  same 
lift  as  the  multiplane  considered).  In  addition,  the  ra- 
tio, K , between  the  induced  drag  of  a multiplane  and 
that  of  the  equivalent  monoplane,  decreases  as  the 
number  of  wings  of  the  multiplane  increases  (with 
global  height  and  lift  kept  constant).  In  the  limit 
as  the  number  of  wings  tends  to  infinity  (“infinity- 
plane”,  according  to  the  Prandtl  [25]  definition),  K 
has  the  minimum  value.  Finally,  there  exists  a box- 
wing type  configuration  that  has  the  same  distribu- 
tion of  circulation  (and  hence  the  same  A")  as  the 
“infinity-plane”;  such  configuration  is  therefore  des- 
ignated, still  by  Prandtl  [25],  at  the  ‘Best  Wing  Sys- 
tem’ (BWS). 

After  the  work  of  Prandtl  [25],  valid  only  for  the  lim- 
ited case  of  multiplanes  with  elliptic  distribution  of 
circulation,  several  numerical  and  experimental  stud- 
ies have  been  published.  Of  particular  interest  here  is, 
the  work  of  Nenadovitch  [23]  on  the  efficiency  of  bidi- 
mensional  biplanes.  Also,  an  experimental  analysis  of 
joined  wing  aircraft  is  presented  by  Wolkovitch  [32], 
whereas  experimental  as  well  as  numerical/empirical 
studies  on  biplanes  and  box-wing  configurations  are 
proposed  by  Gall  and  Smith  [6]. 

In  this  paper,  the  attention  is  focused  on  an  innova- 
tive configuration  for  a new  large  dimension  aircraft 
(usually  indicated  as  New  Large  Airplanes  or  NLA) 
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proposed  by  Frediani  et  al.  [4],  [5].  This  configura- 
tion is  a biplane,  with  counter-swept  wings  (positive 
sweep  for  the  front  lower  wing  and  negative  sweep 
for  the  back  upper  wing,  which  acts  as  a horizon- 
tal tail  as  well)  connected  to  each  other  by  aerody- 
namic surfaces  (Figures  1 and  2).  Following  Frediani 
et  al.  [4],  [5]  we  will  refer  to  this  configuration  as 
the  Prandtl-plane.  Preliminary  numerical  and  exper- 
imental studies  (Frediani  et  al.  [4],  [5])  have  shown 
that  the  induced  drag  of  this  configuration  is  lower 
than  the  induced  drag  of  an  equivalent  monoplane. 
This  fact  allows  one  to  reduce  the  wingspan  of  this 
configuration  without  drag  penalties  and  introduces 
the  possibility  of  respecting  the  maximum  span-wise 
dimensions  (critical  for  the  NLA  with  classical  wing 
configuration),  which  would  allow  compatibility  with 
the  existing  airport  regulations. 

Still  to  in  order  to  put  the  work  in  the  paper  context, 
note  that  through  an  unrelated  activity,  the  authors 
have  been  involved  in  the  developement  of  an  MDO 
code.  The  weakest  module  in  this  code  is  the  aerody- 
namic one.  In  fact,  the  aerodynamic  module  in  such 
a code  is  limited  to  potential  flows  around  simply- 
connected  domains;  on  the  contrary,  the  Prandtl- 
plane  configuration  is  multiply-connected.  In  addi- 
tion, a key  for  this  configurations  (and  MDO  in  gen- 
eral) is  the  drag,  induced  and  viscous. 

These  issues  are  addressed  in  this  paper,  as  a first 
step  towards  generating  a suitable  viscous  aerody- 
namic module  around  multiply-connected  configu- 
rations, to  be  incorporated  in  the  above  mentioned 
MDO  code.  Specifically,  because  of  the  complexity 
and  novelty  of  the  configurations  discussed  above, 
it  is  desirable  to  develop  an  aerodynamic  methodol- 
ogy for  multiply-connected  configurations  capable  of 
yielding  accurate  predictions  with  a relatively  small 
computational  effort  (hence,  not  a CFD  approach), 
which  takes  into  account  the  effects  of  viscosity  and 
of  the  wake  roll-up  (the  geometry  of  the  rolled-up 
wake  is  evaluated  through  the  free-wake  analysis  of 


Ref.  [29]).  The  relatively  recent  state  of  our  work  in 
this  context  was  presented  in  Bernardini  et  al.  [1], 
which  is  based  on  the  classical  viscous-potential  in- 
teraction, with  the  viscous  flow  (in  boundary  layer 
and  wake)  solved  by  using  the  strip-theory  approach, 
with  an  integral  two-dimensional  boundary  layer;  the 
viscous/inviscid  coupling  technique  is  based  upon  the 
Light  hill  [12]  equivalent-sources  approach.  Match- 
ing of  the  boundary-layer  solution  with  the  corrected 
potential-flow  solution  is  obtained  by  direct  iteration. 
Here,  the  emphasis  is  on  the  most  recent  improve- 
ments with  respect  to  that  paper.  First,  we  present 
an  extension  to  higher  order  of  the  direct  boundary 
element  formulation  used  in  the  past  by  the  authors 
for  incompressible  potential  flows  around  lifting  ob- 
jects of  arbitrary  shape  in  uniform  translation;  this 
provides  also  an  opportunity  to  clarify  certain  the- 
oretical issues  connected  with  the  trailing  edge  con- 
ditions. Second,  we  discuss  a three-dimensional  in- 
tegral boundary-layer  formulation,  which  is  based 
essentially  on  the  algorithm  of  Nishida  [24]  and 
Milewski  [15].  The  solution  is  obtained  through  the 
simultaneous  coupling  technique  introduced  by  Drela 
[3], 

2.  Formulation  of  Ref.  [1] 

The  formulation  of  Bernardini  et  al.  [1]  is  outlined 
here  in  order  to  emphasize  the  differences  of  the  new 
formulation  and  also  because,  some  new  results  pre- 
sented here  are  based  on  this  formulation.  For  a de- 
tailed discussion,  see  Ref.  [18],  where  a review  is  also 
presented. 

Potential-Flow  Formulation:  an  inviscid,  incompress- 
ible, initially-irrotational  flow  remains  at  all  times 
quasi-potential  (z.e.,  potential  everywhere  except  for 
the  wake  surface,  which  is  the  locus  of  the  points  em- 
anating from  the  trailing-edge).  In  this  case,  the  ve- 
locity field,  v,  may  be  expressed  as  v = (where 
< p is  the  velocity  potential).  Combining  with  the  con- 
tinuity equation  for  incompressible  flows,  V • v = 0, 
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yields 

vV  = o.  (i) 

The  boundary  conditions  for  this  equation  are  as  fol- 
lows. The  surface  of  the  body,  SB , is  assumed  to  be 
impermeable;  this  yields  (v  — vB)  • n = 0,  i. e. , 

= vB  n for  X G SB,  (2) 

where  vB  is  the  velocity  of  a point  xG<SB)  whereas 
n is  the  outward  unit  normal  to  SB . At  infinity,  in  a 
frame  of  reference  fixed  with  the  unperturbed  fluid, 
we  have  (p  = 0.  The  boundary  condition  on  the  wake 
surface,  Sw\  are:  (i)  the  wake  surface  is  impermeable, 
and  (ii)  the  pressure,  p,  is  continuous  across  it.  These 
imply  that  A (d<p/dn)  = 0 (where  A denotes  discon- 
tinuity across  5^),  and  that  A (p  remains  constant  in 
time  following  a wake  point  xw  (whose  velocity  is  the 
average  of  the  fluid  velocity  on  the  two  sides  of  the 
wake),  and  equals  the  value  it  had  when  xw  left  the 
trailing  edge.  This  value  is  obtained  by  imposing  the 
trailing-edge  condition  that,  at  the  trailing  edge,  A ip 
on  the  wake  equals  pu  — <pi  on  the  body  (subscripts 
u and  / denote,  respectively,  upper  and  lower  sides  of 
the  body  surface). 

In  this  paper,  the  above  problem  for  the  velocity  po- 
tential is  solved  by  a boundary  element  formulation. 
The  boundary  integral  representation  for  the  above 
problem  is  given  by 

'AA  = (ax  - xx~  j dS(y)  (3) 

- I 

w 

with  x :=  VB  • n,  and  G — — l/47r||y  — x||.  Note  that, 
in  the  absence  of  the  wake,  Eq.  3,  in  the  limit  as  x 
tends  to  SB , represents  a boundary  integral  equation 
for  <p  on  SB , with  x on  $b  known  from  the  bound- 
ary condition.  Once  ip  on  the  body  is  known,  p (and 
hence  v and,  by  using  Bernoulli’s  theorem,  p)  may 
be  evaluated  everywhere  in  the  field.  The  situation  is 
similar  in  the  presence  of  the  wake,  since,  by  apply- 
ing the  wake  and  trailing-edge  conditions,  Ap  on  the 


wake  may  be  expressed  in  terms  of  (p  over  the  body 
at  preceding  times. 

In  Bernardini  et  al  [1],  the  boundary  integral  equa- 
tion for  (p  based  on  Eq.  3 is  solved  numerically  by 
discretizing  the  body  and  wake  surfaces  in  quadrilat- 
eral elements,  assuming  py  x and  Ap  to  be  piecewise 
constant,  and  imposing  that  the  equation  be  satis- 
fied at  the  center  of  each  body  element  (zeroth-order 
boundary-element  collocation  method). 

It  should  be  noted  that  the  geometry  of  the  wake  is 
not  known  a priori . For  the  case  of  isolated  wings 
in  uniform  translation,  a flat  wake  with  vortical  lines 
parallel  to  the  unperturbed  streamlines  is  typically 
used  (prescribed-wake  analysis).  In  the  configura- 
tion of  interest  here,  such  a priori  assumptions  are 
not  justified  by  past  experience.  Still  in  Bernardini 
et  al.  [1],  the  wake  geometry  is  either  prescribed 
(undisturbed)  or  obtained  through  a free-wake  anal- 
ysis, z.e the  boundary  integral  equation  is  used  to 
compute,  at  each  time  step  tn  the  flow  velocity  at 
nodes  xt  of  the  discretized  wake  surface  (free-wake 
analysis);  the  locations  of  x,  at  tn+1  is  then  obtained 
by  integrating  the  equation  xt  = v^xj  by  the  ex- 
plicit Euler  method  (in  other  words,  in  the  case  of 
free-wake  analysis  the  shape  of  the  wake  and  the  flow- 
field  solution  are  obtained  step-by-step,  as  part  of 
the  solution).  For  the  steady-state  cases  of  interest 
here,  the  solution  is  obtained  in  the  limit  as  the  solu- 
tion for  t —¥  oo  of  a transient  flow  due  to  an  impul- 
sive start.  Incidentally,  this  procedure  resolves  also 
the  non-uniquess  issue  connected  with  steady-state 
flows  around  multiply-connected  region;  indeed,  in 
the  case  of  unsteady  flow  the  solution  is  unique,  be- 
cause Kelvin’s  theorem  allows  one  to  obtain  the  cir- 
culation on  the  wing  from  the  vorticity  shed  in  the 
wake  (see  Ref.  [20]). 

Boundary  Layer:  the  viscous  flow  analysis  is  limited 
to  attached  steady  high- Reynolds  number  flows.  Un- 
der these  assumptions,  a classical  integral  boundary- 
layer  formulation  may  be  used.  Specifically,  Bernar- 
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dini  et  al.  [1]  use  a 2D  integral  boundary-layer  formu- 
lation used  as  'strip-theory’.  The  laminar  portion  is 
computed  by  using  the  Thwaites  collocation  method 
[30].  The  transition  from  laminar  to  turbulent  flow  is 
detected  by  the  Michel  method  [14].  The  turbulent 
portion  of  the  boundary  layer  and  wake  are  studied 
by  the  dag-entrainment’  method  of  Green  et  al.  [9]. 

Viscous/Inviscid  Coupling:  once  the  boundary-layer 
equations  are  solved,  the  viscous  correction  to  the 
potential  flow  is  obtained  as  a transpiration  velocity 
through  Sb  and  Sw  (equivalent  source  method  by 
Lighthill  [12]).  The  boundary  condition  for  (p  over 
the  body  surface  is  modified  as  follows 

d<p/dn  = vB  - n + Xv. 

where  the  transpiration  velocity  Xv  is  given  by 


where  s denotes  the  arclength  in  the  2D  boundary- 
layer  and  wake,  and  5*  is  the  displacement  thickness. 
In  addition,  u and  ue  denote,  respectively,  the  veloc- 
ity within  the  vortical  layer  and  at  its  outer  edge  (in  a 
frame  of  reference  fixed  with  the  body).  On  the  wake 
surface  one  has 

A (d<p/dn)  = (Xv)u  + (Xv)<» 

where  the  subscripts  u and  l denote,  respectively,  the 
upper  and  lower  sides  of  the  wake  surface. 

3.  Third-order  extension 

The  specific  formulation  used  here  is  based  upon  a 
high-order  formulation  introduced  by  Gennaretti  et 
al.  [8],  and  extended  in  Refs.  [20]  and  [21],  to  which 
the  reader  is  referred  for  details  (here  the  emphasis 
is  on  the  applications).  This  consists  of  using  bicubic 
quadrilateral  elements  and  the  Hermite  interpolation 
for  (p  in  both  directions.  This  yields  a bicubic  repre- 
sentation in  terms  of  the  nodal  values  of  <£>,  d<p/dl;> 
dp> /dr] , and  d2<p/dl;dri\  the  nodal  derivatives  are  then 
expressed  in  terms  of  the  unknowns  of  the  problem, 


(fii,  through  suitable  finite-difference  approximations. 
Only  steady  flows  are  considered  here;  then,  the  wake 
geometry  at  the  trailing-edge  (typically  tangents  to 
one  of  the  sides  of  the  body  surface  at  the  trailing- 
edge)  may  be  determined  by  the  results  of  Mangier 
and  Smith  [13]. 

The  above  formulation  generates  new  issues  con- 
nected to  the  trailing-edge.  The  first  one,  is  related 
to  the  existence  of  the  trailing-edge  potential  discon- 
tinuity which  implies  that  at  any  trailing-edge  node 
there  exist  two  unknowns,  but  only  one  collocation 
point.  Gennaretti  et  al.  [8]  use  two  collocation  points 
suitably  located  slightly  ahead  of  the  trailing  edge 
(otherwise  the  collocation  points  coincide  with  the 
nodes). 

In  order  to  avoid  ill-conditioned  matrix  here,  as  men- 
tioned above,  we  use  the  approach  introduced  in  Ref. 
[20],  who,  at  the  trailing-edge,  use  only  one  control 
point  but  two  different  integral  equations:  the  first  is 
the  potential  equation  while  the  second  is  its  deriva- 
tive in  the  direction  of  the  wake  normal  at  trailing- 
edge. 

An  additional  issue  stemming  from  the  third-order 
formulation  is  that  of  the  evaluation  of  the  derivatives 
of  <p  at  the  trailing  edge.  These  may  be  obtained  by 
imposing  the  two  conditions  of  stagnation  point  and 
smooth  flow  (see  Ref.  [20]  for  details). 

4.  Three-dimensional  boundary  layer 
In  order  to  take  into  account  the  three-dimensional 
effects,  the  viscous  flow  in  the  boundary  layer  and 
wake,  is  modelled  using  the  three-dimensional  in- 
tegral boundary-layer  equations  on  non-ortogonal 
grids.  Specifically,  the  integral  method  uses  two  equa- 
tions for  momentum  (extention  of  von  Karman  equa- 
tion to  3D  flows)  coupled  by  two  auxiliary  equations. 
The  first  is  the  kinetic  energy  equation  and  the  sec- 
ond is  the  transport  equation  for  the  maximum  shear 
stess  coefficient  ('lag’  equation).  In  order  to  com- 
plete the  problem,  we  use,  following  Nishida  [24]  and 
Milewski  [15],  the  streamwise  closure  relations  pro- 
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posed  by  Drela  [3]  in  his  2D  integral  boundary-layer 
method,  along  with  the  crosswise  relations  of  John- 
ston [11].  The  Mughal  [22]  finite-volume  scheme  is 
used  for  the  solution. 

Also  in  this  case,  the  viscosity  correction  to  the  po- 
tential flow  is  evaluated  as  a transpiration  velocity 
through  Sb  and  Sw  (equivalent  source  method  by 
Lighthill  [12]).  However,  as  mentioned  above,  in 
contrast  with  the  formulation  of  Bernardini  et  ol. 
[1],  the  matching  of  the  thin-layer  solution  with  the 
potential-flow  solution,  is  obtained  by  the  fully  simul- 
taneous coupling  method  of  Drela  [3],  which  solves 
the  viscous  and  inviscid  equations  simultaneously  and 
is  stable  even  for  separated  flows. 

5.  Numerical  results  and  comments 

The  results  presented  here  fall  into  three  groups.  In 
the  first  group  are  some  new  results  obtained  with 
the  formulation  outlined  in  Section  2 (some  of  the 
results  of  Ref.  [1]  are  included  here  for  the  sake  of 
clarity).  The  second  and  third  group  are  based  on 
the  formulations  of  Sections  3 and  4. 

For  all  the  potential  flow  results,  the  airloads  are  de- 
termined by  using  the  formulation  of  Ref.  [7],  which 
is  an  exact  generalization  of  the  Trefftz-plane  theory 
[31]  for  the  evaluation  of  aerodynamic  loads  around 
objects  of  arbitrary  shapes.  The  viscous  drag  is  eval- 
uated by  using  the  classical  formula  by  Squire  and 
Young  [28]. 

First,  a few  results  obtained  by  Bernardini  et  al.  [1] 
are  discussed  as  they  help  in  putting  the  paper  in  the 
proper  context.  First  of  all,  the  potential  flow  for- 
mulation of  Section  2 has  been  validated  in  Ref.  [1] 
by  a comparison  with  analytical  results  by  Prandtl 
[25]  either  in  the  case  of  biplanes  and  box- wing  con- 
figuration; this  comparison  has  shown  nearly  perfect 
agreement  between  our  numerical  results  and  the  an- 
alytical results  of  Prandtl.  Also,  in  order  to  analyze 
the  mechanism  of  induced-drag  reduction,  it  is  inter- 
esting to  compare  the  efficiency  K of:  (i)  a rectangu- 


lar biplane;  (it)  a rectangular  biplane  with  facing  tip 
winglets  on  the  two  wings  and  (Hi)  a box-wing.  For, 
this  show  the  graduality  of  the  reduction  (in  the  sense 
that  the  biplane  corresponds  to  the  zero  winglet- 
length  case  and  the  wing-box  to  the  full  length  case); 
this  is  apparent  from  Figure  3 (from  Ref.  [1]),  which 
depicts  a comparison  of  the  curves  K = K (G/b)  for  a 
biplane  without  winglets,  for  the  corresponding  box- 
wing body,  and  for  a biplane  with  winglets  (winglet 
length  0.4c,  where  c is  the  root  chord). 

It  should  be  observed  that  the  above  results  have  been 
obtained  using  a prescribed-wake  approach,  in  or- 
der to  be  consistent  with  the  Prandtl  model  which  is 
based  on  the  lifting-line  theory.  However,  the  present 
methodology  is  able  to  capture  the  effects  of  the  wake 
roll-up  on  the  airloads  by  using  a free-wake  approach, 
as  described  above.  Figure  4 (also  from  Ref.  [1])  de- 
picts the  Prandtl  factor  K as  a function  of  Gjb)  for 
free-wake  and  prescribed-wake  analysis.  The  results 
obtained  by  a free-wake  approach  show  that  the  wake 
roll-up  effect  is  to  reduce  the  value  of  K.  For  the 
range  of  G/b  of  practical  interest  in  aeronautical  ap- 
plications (he.,  0.1  < G/b  < 0.2),  we  have  that  the 
reduction  of  K is  about  6%  for  the  case  of  a box- 
wing. This  demonstrates  the  importance  of  including 
a rolled-up  wake  in  the  analysis. 

Next,  consider  new  potential-flow  results  obtained 
with  the  formulation  of  Ref.  [1].  In  particular,  Fig.  5 
present  a parametric  study  of  the  effects  of  the  sweep 
angle  (of  the  front  and  rear  wing)  on  the  Prandtl  fac- 
tor K of  the  Prandtl-plane  shown  in  Figures  1 and 
2.  Here,  we  have  considered  three  different  config- 
urations characterized  by:  (a)  Mid-chord  front  wing 
sweep  angle  Ap  = 35°  and  mid-chord  rear  wing  sweep 
angle  A r = —40°,  (b)  mid-chord  front  wing  sweep 
angle  Ap  = 37°  and  mid-chord  rear  wing  sweep  an- 
gle A r = —37°  and,  (c)  mid-chord  front  wing  sweep 
angle  Ap  = 40°  and  mid-chord  rear  wing  sweep  an- 
gle A r = —35°.  We  can  note  that  the  Prandtl  factor 
of  all  three  configuration  are  hardly  distinguishable, 
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and  are  very  close  to  the  value  of  the  Prandtl  factor 
of  the  Best  Wing  System  (continuous  line).  This  is 
very  important  because  it  shows  that,  even  for  the 
Prandtl  plane,  we  have  a reduction  of  induced  drag, 
with  respect  to  a traditional  wing  configuration. 
Next,  consider  some  preliminary  results  for  attached 
high-Reynolds  viscous  flows.  The  two-dimensional 
integral  boundary-layer  formulation,  used  a ^trip- 
theory’  in  three-dimensional  application,  has  been 
validated  in  Ref.  [1],  by  comparison  with  experimen- 
tal results  available  in  literature,  in  the  case  of:  (i) 
isolated  wing,  (it)  biplane,  and  (Hi)  box-wing  config- 
uration. For  the  sake  of  completeness  the  results  for 
the  third  case  are  presented  here  in  Fig.  6 (also  from 
Ref.  [1])  which  depicts  the  polar  at  Re  = 5.1  * 105 
of  a box-wing  configuration  with  aspect  ratio  5,  gap 
G = c,  stagger  S = c and  section  profiles  NACA  0012. 
The  numerical  results  are  in  good  agreement,  also  in 
this  case,  with  the  experimental  results  by  Gall  and 
Smith  [6].  Next  consider  new  results.  The  formula- 
tion of  Ref.  [1]  has  been  used  for  a parametric  study 
of  the  effect  of  sweep  angle  on  the  polar  of  a Prandtl- 
plane.  This  is  shown  in  Fig.  7 which  depicts  a po- 
lar for  the  same  three  configuration  of  Prandtl-plane 
used  in  Fig.  5.  We  can  note  that  their  difference  is 
very  small. 

Next  consider  the  validation  of  third-order  formula- 
tion. Some  results  have  been  presented  in  Morino  and 
Bernardini  [21]  for  the  case  of  simply  (RAE  wing)  and 
multiply  connected  (nacelle)  configurations.  Here,  we 
consider  the  particularly  challanging  case  of  a wing 
with  strake  (see  Fig.  8),  at  angle  of  attack  a = 5°  and 
having  symmetric  NACA  0002  airfoil  sections.  Fig. 
9 depicts  the  circulation  F = A <pTE  as  a function 
of  2 y/b)  while  Fig.  10  shows  the  chordwise  distri- 
bution of  the  pressure  coefficient  (at  2 y/b  = .219). 
The  present  results  are  obtained  with  a chordwise 
discretization  N\  ~ 20  and  a spanwise  discretization 
N2  = 12,  and  are  compared  with  those  obtained  with 
the  method  of  Roberts  and  Rundle  [26]  with  N\  = 39 


and  N2  = 24  and  those  obtained  with  the  method  of 
Rubbert  and  Saaris  [27]  with  N\  = 28  and  N2  = 12. 
The  agreement  among  the  three  formulations  is  good, 
except  at  the  trailing-edge,  where  the  pressure  pre- 
dicted in  our  case  is  higher  (see  also  Ref.  [21]  which 
shows  why  our  results  appear  to  be  more  reliable). 
Finally,  we  present  some  preliminary  results  of  the 
three-dimensional  boundary  layer  formulation.  In 
particular,  Figs.  11  and  12  show  respectively  the 
chordwise  distribution  (at  77  = 0.5  and  for  the  upper 
surface  of  the  wing)  of  the  streamwise  displacement 
thickness,  <$*,  and  the  shape  factor,  H , for  a RAE 
wing  with  mid-chord  sweep  angle  A = 45°,  RAE 
101  airfoil  secton  at  Re  — 2.1  ■ 106  and  angle  of  at- 
tack a — 6.3°.  Our  numerical  results,  obtained  using 
Ni  — 60  and  N2  — 10,  have  been  compared  with 
the  numerical  results  of  Milewski  [15]  (obtained  us- 
ing Ni  = 64  and  N2  = 10),  and  experimental  results 
of  Brebner  and  Wyatt  [2].  The  transition  location, 
in  this  case,  is  fixed  using  experimental  data.  We 
have  also  included  in  the  figure  the  two-dimensional 
boundary  layer  formulation  (unfortunately  with  tran- 
sition obtained  by  Michel  method  [14]).  Our  three- 
dimensional  numerical  results  are  in  good  agreement 
with  both  the  Milewski  [15]  and  with  experimental 
data.  On  the  contrary,  the  2D  formulation  shows,  at 
the  trailing  edge,  a great  difference.  This  could  find  a 
justification  in  the  fact  that,  at  the  trailing  edge,  the 
3D  effects  are  higher. 

To  conclude  the  paper  we  wish  to  mention  a tech- 
nique that  allows  us  to  overcome  the  boundary  layer 
limitations.  This  is  based  on  a decomposition  in- 
troduced by  Morino  [17]  for  the  velocity  field  as 
v = where  w (obtained  by  direct  integration 

of  V x w = ( :=  V x v)  has  the  property  that  w = 0 
where  £ = 0.  This  decomposition  has  been  success- 
fully applied  to  2D  steady  full  potential  flows  with 
boundary-layer  correction  (see  Morino  et  al.  [19]) 
and  to  2D  steady  Euler  flows,  obtained  as  an  exact 
correction  to  the  full-potential  solver  (see  lemma  et 
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al  [10]).  The  coupling  of  Euler  with  the  boundary 

layer  is  now  under  investigation. 
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Fig.  1.  Isometric  view  of  Prandtl-plane. 


Fig.  2.  Isometric  view  of  wing  of  Prandtl-plane. 


Fig.  3.  Effect  of  winglets  on  Prandtl  factor:  rectangular 
biplane  (6/c=  10,  a — 5°,  unstaggered). 


G/b 


Fig.  5.  Effect  of  sweep  angle  oil  Prandtl  factor:  Prandtl- 
plane. 


Fig.  4.  Effect  of  wake  roll-up  on  Prandtl  factor:  box-wing 
body  (fc/c  = 10,  a = 5°,  unstaggered). 


Fig.  6.  Polar  of  box-wing  body  (aspect  ratio  5,  G/c  = 1, 
stagger  = c,  section:  NACA  0012)  at  Re  = 5.1  • 105. 
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DISCUSSION 


Session  II,  Paper  #16 


Herr  Sacher  (DASA,  Germany)  asked  how  the  engine  would  be  represented  in  the  potential 
approach. 

Prof  Morino  noted  that  in  the  past,  the  engine  has  been  represented  by  a prescribed  flow 
at  the  inlet  with  a rigid  surface.  As  this  rigid  surface  should  be  coincident  with  the  actual 
stream  surface,  iteration  is  required  to  determine  this  shape.  In  this  present  work,  the 
approach  is  to  model  the  jet  with  a prescribed  uniform  flow  plus  a wake  (doublet  layer)  to 
model  the  velocity  discontinuity  between  the  jet  and  the  surrounding  flow.  It  was  noted 
that  this  would  be  an  ideal  application  of  the  viscous  flow  extension  discussed  during  the 
presentation. 


